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Abstract. The elliptic flow {V2) pattern in terms of hadron mass and transverse 
momentum pT is qualitatively described for pT < 2 GeV/c by ideal hydrodynamics 
in Au + Au collisions at RHIC. In addition, for pT = 2 — 6 GeV/c the 
measured V2 follow a universal scaling by the number of quarks explained by 
quark coalescence/recombination models. These observations suggest that a partonic 
collectivity develops in the matter in early stage of heavy ion collisions. Centrality as 
well as system size and energy dependence of the V2 is important to shed light on the 
underlying collision dynamics in heavy ion collisions. Wc present the measurements of 
centrahty dependence of V2 at ^/s^n = 200 and 62.4 GcV in Au + Au and Cu + Cu 
collisions for Kg, cj), A, S and at STAR experiment. We focus on the recent Cu + 
Cu results and discuss the centrality dependence of V2 as well as the number of quark 
scaling as a function of transverse kinetic energy at different system size and energies. 
We also discuss the eccentricity scaled V2 for identified hadrons and implications that 
ideal hydrodynamical limit has not been reached at RHIC. 
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1. Introduction 

Elliptic flow is expected to be one of the most sensitive observable to study early stage of 
heavy ion collisions at Relativistic Heavy Ion Collider (RHIC), see recent review in [1]. 
The elliptic flow is defined by the second harmonic Fourier coefficient of azimuthal 
distribution of produced particles with respect to the reaction plane 

1^2 = (cos(20-2^Rp)), (1) 

where is the azimuth of particles in the laboratory frame, \&rp is the azimuthal 
angle of reaction plane which is determined by the direction of impact parameter and 
beam axis [2]. Brackets denote the average over all particles and events. In Au + Au 
collisions at RHIC, it has been found that heavier hadrons has smaller V2 than lighter 
hadrons for transverse momentum < 2 GeV/c, which is qualitatively explained by 
hydrodynamical models [3]. Above px = 2 GeV/c, the V2 had a reversed trend where 
protons and A's have larger V2 than the charged mesons [11[5]. The reversed trend can 
be well described by quark recombination/coalescence models [6l[7l[H] that indicate the 
f 2 has already developed during the early partonic stage prior to the hadronization. It 
is important to study the systematics of ^2, such as system size, beam energy as well 
as centrality dependence, in order to understand the underlying collision dynamics in 
heavy ion collisions. 



2. Data Analysis 

In this analysis, we used minimum bias 24 M events in Cu + Cu collisions at 
■^snn = 200 GeV taken at STAR experiment. The event centrahty was determined 
by the uncorrected charged particle multiplicity measured by the main TPC within 
pseudorapidity |?7| < 0.5. Particle identifications were performed by the energy loss in 
the main TPC that enable us to separate Kg, A and H in 0.2 < < 4 GeV/c. Event 
plane method [2] was used to measure V2 at midrapidity. The V2 can be expressed as 

^ (cos(20-2v^n) 
(cos(2vl/fP-vl;j,p))' 

where \l/f^ denote the second harmonic event plane determined from elliptic flow, and 
brackets represent the average over detected particles and events. Event plane was 
determined by using the forward time projection chamber (FTPC) covered at 2.5 < \ri\ < 
4. Since the FTPC has a rapidity gap from the main TPC region, it can reduce the 
non-flow contributions that are not correlated with the reaction plane, such as di-jets 
and resonance decays etc. The denominator in ([2]) is defined as an event plane resolution 
that is about 0.18 in midcentral Cu + Cu collisions. We also analyzed the data in Cu 
+ Cu collisions at ^/sNN = 62.4 GeV, which has 12.5 M minimum bias events. The 
event plane resolution was by a factor of 2 smaller than that at 200 GeV due to lower 
multiplicity. 
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3. Results 
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Figure 1. (a) The V2 as a function of for K%, A-H A, S-hS in - 60 % (top), - 20 
% and 20 - 60 % centralities (bottom). Dashed lines represent ideal hydrodynamical 
calculation [9]. Only statistical errors arc shown, (b) Number of quark scaling of V2 
as a function of mx — mass in - 60 % (top), - 20 % and 20 - 60 % (bottom). Data 
symbols are the same as those in (a). 



Figure [H (a) shows the f 2 for K^, A -|- A and S -|- S as a function of in different 
centrahty selections. We have observed that A has smaller V2 than below = 2 
GeV/c in Cu + Cu collisions. For px > 2 CeV/c, however, A V2 become larger than 
that of K^. We have also found S has sizable V2 in minimum bias - 60 % centrality. 
Results were compared to an ideal hydrodynamical calculation with the first order phase 
transition from the QGP to the hadron phase. The ideal hydrodynamical model does not 
describe the centrality dependence of our data. For - 20 %, the model under-predicts 
our data and for 20 - 60 %, it over-predicts the V2. Effects not included in the model 
which may be relevant are geometric fluctuations in the initial conditions (particularly 
important in central collisions) and finite viscosity effects. It remains to be seen if these 
effects can account for the difference between the models and data. Figure [1] (b) shows 
the number of quark (NQ) scaling of V2 as a function of transverse mass minus particle 
mass mx — mass, namely both vertical and horizontal axes are divided by the number 
of valence quarks for each hadron. One can see that all measured hadrons lie on the 
universal curve up to px = 4 GeV/c in all centrality bins. This results indicate that the 
V2 has already established in the partonic stage in Cu + Cu collisions. 

Figure [2] show the comparison of V2 for Kg and A -|- A at y^s^N = 200 and 62.4 
CeV in Cu + Cu collisions (left panel). For comparison, results in Au + Au collisions 
are presented in right panel [101 IH] • We found that the measured V2 at ^/sNN = 200 
GeV are consistent with those at yJwTV = 62.4 GeV in both Cu + Cu and Au + Au 
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Figure 2. Comparison of 112 (pt) for and A + A in - 60 % ccntrality at ^snn 
= 200 GcV with those at 62.4 GeV in Cu + Cu (left) and Au + Au colhsions (right). 
The V2 in Au + Au coUisions are taken from [101 [TT]. Only statistical errors are shown. 
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Figure 3. Number of quark 
scaling of V2 divided by participant 
eccentricity Spart for K'g and A + A 
0jS cL function of mx — mass in Cu 
+ Cu collisions. For comparison, 
results in Au + Au collisions are 
also plotted. 
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collisions. 

Figure [3] shows the eccentricity scaled V2 for and A + A as a function of 
mx — mass. In order to compare different particle species, NQ scaling was also applied. 
Since the event plane at the FTPC was determined by particles from participant 
nucleons, appropriate initial geometrical anisotropy should be participant eccentricity 
Spart that include event-by-event position fluctuation of participants |T2]. The scaled 
V2 in - 20 % appeared to be larger than those in 20 - 60 %, which indicate that 
stronger collective flow has developed in more central collisions. In addition, scaled 
V2 in Cu + Cu collisions are consistent with that in Au + Au collisions with similar 
number of participant Npart- This result suggest the observed V2 is determined not only 
by initial geometry, but also by the Npart- Therefore, it is important to understand 
whether the V2 in most central Au + Au collisions reach ideal hydrodynamical limit or 
not. Because ideal hydrodynamical models predict V2 is only determined by the initial 
spatial anisotropy. 
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Figure 4. (Left) V2/S as a function of transverse number density {\/S)dN/dy for 
different particle species from STAR and PHENIX. STAR results include preliminary 
f2{4}, W2{ZDCSMD} in Au + Au collisions as well as V2 for ifg and A in Cu + Cu 
collisions, and published V2{2} [H], V2 for and A g]. PHENIX results were taken 
from [14]. Error bars denote quadratic sum of statistical and systematic uncertainties. 
(Right) V2/e divided by extracted hydrodynamical limit from the fit as a function of 
transverse number density {\ / S)dN / dy . Dashed line represent the simultaneous fit for 
all measured V2 by ([3]). 



Figure m show the comparison of V2/S as a function of transverse number density 
{1/ S)dN / dy for different particle species in Au + Au colhsions at y/s^N = 200 GeV 
from PHENIX and STAR experiments. Published STAR and PHENIX results were 
taken from [H [13], [14] . In order to extract the ideal hydrodynamical limit, we have used 
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the formula based on a transport model with a finite Knudsen number [15] 

e~\e /Hydro VI + K/Ko) ' K ~ S dy ' ^ ^ 

where ('y2/£)Hydro ideal hydrodynamical limit of V2, K is the Knudsen number, 

Cs is the speed of sound, a is the partonic cross section, S is the transverse area of 
collision zone and dN/dy is the total hadron multiplicity. We adopted Kq = 0.7 and 
Cs = 1/ VS in order to reproduce the transport model calculation [16] . Hydrodynamical 
limits were extracted for different hadrons by fitting the available V2 simultaneously, 
where we assumed a is the same for all hadrons. Right panel in Figure H] shows the ratio 
of V2 to the resulting hydrodynamical limit from the fit as a function of {1/ S)dN/dy. 
The ratio should be 1 when V2 reach hydrodynamical limit. We found that the ratio 
converge into a universal curve and ideal hydrodynamical limit has not been reached 
even at the most central collisions at RHIC. 



4. Conclusions 

In summary, we have measured f 2 for Kg, A + A and S + S at y^Sivjv = 200 GeV in Cu 
+ Cu collisions. We found that A has smaller V2 than Kg below = 2 GeV/c, whereas 
the trend was reversed above = 2 GeV/c. The measured hadrons follow the number 
of quark scaling as a function of mx — mass in Cu + Cu collisions at ^/sNN = 200 GeV. 
This result is consistent with the idea that collective fiow has been established in the 
partonic stage before the hadronization takes place. The V2 at ^/JNN = 62.4 GeV were 
consistent with those in ^/sNN = 200 GeV in Cu + Cu collisions. We observed that 
eccentricity scaled V2 was larger in more central events in Cu + Cu collisions, which 
indicate that stronger collective fiow was developed in more central collisions. Finally, 
from the comparison of the data with a transport model approach, we found that ideal 
hydrodynamical limit has not been reached even at the most central collisions at RHIC. 
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